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Abstract: Employing the recently developed 
NEXus model, we compare the yields of differ- 
ent hadrons in ultra-relativistic collisions: electron- 
positron (e + e~) annihilation at y/s —91 GeV, proton- 
proton (pp) scattering at y/s =17 GeV and nucleus- 
nucleus (AA) collisions at y/s =17 GeV (SPS) and 
yfs =200 GeV (RHIC). Plotting the yields as a func- 
tion of the hadron masses, we find very surprising 
results: we observe that the spectra are practically 
identical for e + e~at 91 GeV and central nucleus- 
nucleus reactions at SPS and RHIC energies, whereas 
the spectrum for proton-proton scattering is some- 
what steeper. All have the form one expects if the 
particles were emitted by a canonical system which is 
characterized by a temperature and chemical poten- 
tials. These identical forms have, however, different 
origins: in e + e~and pp the exponential shape it is 
due to the statistical behavior of string fragmenta- 
tion, which has absolutely nothing to do with ther- 
malization, in AA it is caused by phase space. The 
fact that e + e~and nuclear results agree is pure coinci- 
dence. Surprisingly the results for pp and e + e~differ, 
although here the production mechanism is identical. 
In pp collisions we see directly that the string energy 



is very limited and hence the high mass baryons are 
suppressed. We conclude that it is practically impos- 
sible to draw conclusions from hadronic yields about 
the reaction mechanism. 



1 Introduction 

Studying hadron production in e + e~ annihilation and 
pp scattering has already a long traditionjlj In 
recent times particle production has gained renewed 
interest because it has been argued that in heavy 
ion collisions the enhanced yield of strange particles 
(as compared to pp) may serve as "smoking gun" 
for the conjectured creation of a quark gluon plasma 
(QGP)Q]. Looking carefully at the presently avail- 
able data of the reaction 158 AGeV Pb + Pb one 
observes indeed this fact Q . 

The formation of particles with a given mass de- 
pends always on the energy which is available in the 
system which forms them. If this energy is not much 
higher than the masses of the final state particles one 
expects that the phase space of the strange particles 
is much smaller than the one of their lighter non- 



strange counterparts and therefore they are less fre- 
quently produced. This is expected to be the case in 
two body collisions in an expanding hadron gas with 
a temperature of T=170 MeV as considered as a good 
description of the final state of the heavy ion inter- 
action. Contrary, if the energy of the forming system 
is much larger than the hadron masses, the mass dif- 
ference should not play an important role and we ex- 
pect to see (almost) the same yield. This scenario is 
expected if the hadrons are produced by a decaying 
QGP. 

In this situation it seems to be worthwhile to apply 
state-of-the -art models in order to check the validity 
of the above mention conjecture. Here we use the 
recently developed NEXUS model, whose results have 
been extensively compared with most of the existing 
e + e~, pp, fj,A, pA and AA data 0], |f|. For our study 
we compare the results of simulations for 4 different 
reactions: e + e~at ^=91 GeV, pp at y / s=17 GeV, 
central Pb-Pb collisions at \/s=17 AGeV (SPS) and 
central Au-Au collisions at ^5=200 AGeV (RHIC). 

2 NEXUS 

How can one simulate in a realistic way hadron- 
hadron or even nucleus-nucleus collisions? A sophis- 
ticated approach to such high energy interactions is 
the so-called Gribov-Regge theory. This is an ef- 
fective theory, which allows multiple interactions to 
happen "in parallel", with phenomenological object 
called "Pomerons" representing elementary interac- 
tions. Using the general rules of field theory, on may 
express cross sections in terms of a couple of param- 
eters characterizing the Pomeron. Interference terms 
are crucial, they assure the unitarity of the theory. 
A big disadvantage of this approach is the fact that 
cross sections and particle production are not calcu- 
lated consistently: the fact that energy needs to be 
shared between many Pomerons in case of multiple 
scattering is well taken into account when calculat- 
ing particle production (in particular in Monte Carlo 
applications), but not for cross section calculations. 
Another problem is the fact that at high energies one 
needs a consistent approach to include both soft and 
hard processes. The latter ones are usually treated 



in the parton model approach, which only allows to 
calculate inclusive cross sections. 

We recently presented a completely new approach 
|j| for hadronic interactions and the initial stage of 
nuclear collisions, which is able to solve several of 
the above-mentioned problems. We provide a rigor- 
ous treatment of the multiple scattering aspect, such 
that questions as energy conservation are clearly de- 
termined by the rules of field theory, both for cross 
section and particle production calculations. In both 
(!) cases, energy is properly shared between the 
different interactions happening in parallel. This 
is the most important and new aspect of our ap- 
proach, which we consider a first necessary step to 
construct a consistent model for high energy nuclear 
scattering. The elementary interaction is the sum of 
the usual soft Pomeron and the so-called semi-hard 
Pomeron, where the latter one may be obtained from 
perturbative QCD calculations (parton ladders). To 
some extend, our approach provides a link between 
the Gribov-Regge approach and the parton model, 
we call it "Parton-based Gribov-Regge Theory". Al- 
though the main purpose of NEXus is the simulation 
of nuclear collisions, also electron-positron annihila- 
tion and deep inelastic lepton-nucleon scattering have 
been intensively studied to make sure to have a con- 
sistent approach of high energy scattering. 

The NEXus simulation program is therefore 
available to simulate lepton-lepton, lepton-hadron, 
hadron-hadron, hadron-nucleus, and finally nucleus- 
nucleus collision at high beam energy. 

In electron-positron annihilation, starting from the 
virtual Z, a time-like cascade of parton develops, 
which are used to define a kinky string via identi- 
fying partons and kinks. Particle production is the 
calculated using the standard covariant string frag- 
mentation procedures: the string breaks into string 
segments vis quark-pair production, the segment rep- 
resenting the final hadrons and resonances. 

In pp scattering, one has to deal with a compli- 
cated multiple scattering process, but at the end the 
situation is quite similar to electron-positron anni- 
hilation in the sense that parton configurations are 
translated into string language in order to "parame- 
terize" the theoretically unknown hadronization. 

Concerning nucleus-nucleus collisions, we have to 



2 



distinguish primary and secondary interactions. Pri- 
mary interactions are a straightforward generaliza- 
tion of proton-proton scattering: multiple scattering 
- parton production - string formation. Whereas in 
pp only few strings are formed and one may with- 
out problem consider independent string fragmenta- 
tion, this is not realistic for nuclear collisions. This 
fact is take into account by so-called secondary in- 
teractions, going beyond neXus. To get some idea 
what happens we employ a simple droplet model: 
based on neXus, we calculate energy densities at a 
given proper time (say at 1 fm/c). Spatial region with 
an energy density bigger than some critical value are 
considered to "quark matter" and we treat the further 
development macroscopically, assuming Bjorken-like 
longitudinal expansion and some transverse expan- 
sion. Finally, when the energy density drops below 
the critical density, the droplet decays into hadrons, 
purely via phase space. This approach proved to 
provide a very satisfactory description of the CERN 
heavy ion reactions M. 
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3 Results 



Before embarking to detailed comparisons and con- 
clusions, we have to make sure that the model de- 
scribes reasonable well the observables we use here. 
Here we are interested in the multiplicity in 4ir and 
therefore it is useful to compare the Nexus results 
with experiments as well in 4tt J?], [|, In fig. [l], we 
compare the 4ir Nexus results for particle multiplic- 
ities (after decay of all not identified resonances) for 
central events with a compilation of the experimental 
data. The Q yield was obtained by assuming that the 
ratio of S to f2 at mid-rapidity obtained by WA97 is 
the same as the one for Air. On the right hand side, 
we display mesons and baryons, on the left hand side 
mesons and anti-baryons. The Nexus results are in 
good agreement with the extrapolated experimental 
data over four orders of magnitude. In view of a sim- 
ilar good agreement for e + e~and pp || it is useful 
to use Nexus for the question at hand. In fig. |l|, we 
display at the bottom our predictions for central col- 
lisions of Au + Au at y/s = 200 AGeV (here as well 
non identified resonances have decayed). 
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Figure 1: Top row: The multiplicity per event of 
selected hadrons in 4ir as a function of their mass 
as compared to the experimental values which have 
been extrapolated to 4tt. Bottom row: Predicted 
multiplicity (after resonance decay for the RHIC ex- 
periment Au+Au at GeV. On the left hand side we 
display mesons and baryons, on the right hand side 
mesons and anti-baryons. 



3 



pp 17 GeV 




Figure 2: Distribution of the string energies in 
e+e"at ^fs = 9lGeV , pp and Pb+Pb at y/s = 17 
GeV and Au+Au at Vs=200 GeV. 

The important difference of pp and AA as com- 
pared to e + e~reactions is the distribution of the 
energy W of the strings, i.e. the energy which is 
available for the production of particles. For the 
4 different reactions it is displayed in fig.||. In the 
e + e~reactions at 91 GeV gluon bremsstrahlung rarely 
produces a quark-anti-quark pair and hence the de- 
caying string has almost always the full energy of 91 
GeV. In pp and AA reactions, on the contrary, only 
a small fraction of the initial energy of the protons 
is available in the strings. The distribution of the 
string energies for the reactions at *Js = 17 GeV has 
a peak at around lGeV and falls steeply off towards 
the maximal value of 17 GeV. The string contains 
only a small fraction W = y / xiX2S of the available 
energy -y/s, due to the fact that the light cone mo- 
mentum fractions x\ and x-i of the partons represent- 
ing the string ends are effectively given by structure 
functions which peak at low x. This peak value is 
small as compared to the threshold for multi strange 
baryon production (due to strangeness conservation 
the threshold for ft is around 3 GeV). For RHIC en- 
ergies the average string energy is already very large 
as compared to the mass of stable baryons but still 
smaller than that for e+e~. 

Due to these distributions of string energies, we ex- 
pect to see differences in the yield of heavy as com- 
pared to light hadrons in the different reactions. In 
figj^, we display the yields for the different reactions 
divided by m( 3 / 2 )(2J +1). In the case of AA col- 
lisions, only those detected particles are considered 



which have been created directly from string decay 
and not via droplets. In any case, we observe a 
roughly exponential drop of the yields as a function of 
the mass: it is always more difficult to produce high 
mass particles compared to low mass ones. The same 
exponential form of the mass yield is obtained in this 
representation in canonical models which describe the 
mass yield in terms of a temperature and chemical po- 
tentials |fLol| . Looking more closely we observe, how- 
ever, different slopes of theses exponential functions: 
pp (17 GeV) and A A (17 GeV) show a similar behav- 
ior and both drop faster than e + e~(91 GeV) and AA 
(200 GeV). In order to see the effects more clearly, 
we divide all distributions by the e + e~results, see 
fig. |. The e + e"(91 GeV) - by definition - and AA 
(200 GeV) show a horizontal spectrum, whereas pp 
(17 GeV) and A A (17 GeV) show both a suppression 
of heavy particles as compared to e + e~ . At 17 GeV 
the string masses are small and hence phase space re- 
duces the production of heavy particles, whereas for 
e + e"(91 GeV) and AA (200 GeV) there is no such 
limitation. 

So far we considered only AA collisions with 
hadrons coming directly from strings. The majority 
of hadrons which are observed finally is, however, pro- 
duced from droplets. These droplets decay according 
to the microcanonical phase space into baryons (the 
spin 3/2 decuplet and the spin 1/2 octet) and mesons 
(pseudoscalar and vector octet) respecting all the 
conserved quantum numbers of the droplet ||. Their 
energy is usually large as compared to the hadron 
masses. In fig. we combine the hadrons produced 
by droplets and string decays and show the full AA 
results as seen by the detectors, again compared to pp 
and e + e~ , and in all cases divided by the e + e~results. 
We observe that all the nuclear results (AA (17 GeV) 
and AA (200 GeV)) are almost parallel, whereas the 
pp spectrum drops much faster. If we take pp as 
a reference, we might conclude: we observe indeed 
an enhancement of heavy (multi-strange) hadrons for 
nuclear collisions but also for e + e~, however, in view 
of the above discussion it is more appropriate to talk 
about a strangeness suppression due to the limited 
energy in pp. 

This is a very astonishing result in view of the fact 
that the origin of the hadron production is very differ- 
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Figure 3: Multiplicity /event in 4ir as a function of the 
hadron mass for the 4 investigated reactions. On the 
left hand side we display the mesons and baryons, on 
the right hand side the mesons and anti-baryons. The 
multiplicity has been divided by (2J+1) and m 3 / 2 
(see text). Feeding from higher lying resonances has 
been excluded. In case of AA collisions, we consider 
only those hadrons which are created directly from 
string decay (not via droplets). 



Figure 4: Same as fig ||, but divided by the multiplic- 
ity observed in e + e~at 91 GeV. 



Figure 5: Same as fig but showing the full result for 
AA collisions, i.e. hadrons coming either from string 
decay or from droplets. 
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Figure 6: Multiplicity /event in 4ir as a function of the 
hadron mass for e + e~and droplet of 10 GeV assuming 
different freeze out energy densities e. 



ent in e + e~as compared to the hadronic interactions. 
As said before, the decay of a string in e + e~is deter- 
mined by the area law which is close to a sequential 
decay in a one dimensional longitudinal phase space. 
The droplet, on the contrary, disintegrates instan- 
taneously and the different decay channels are popu- 
lated according to the three dimensional phase space. 
Is it a accident that both mechanism give about the 
same hadron multiplicity distribution? The answer 
is clearly yes as can be seen from fig. ^. 

The droplet decay pattern depends on the energy 
density e where the decay takes place. Smaller energy 
density implies larger volumes and therefore decay 
channels in which many particles are produced are 
favored due to the additional volume factor for each 
additional particle. If more particles are produced 
the heavier ones are more suppressed, as can be seen. 
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Now it happens to happen that at energy densities 
of e^lGeV, a number which can be derived on the- 
oretical grounds as well as by comparison with the 
data, e + e~and hadronic interactions have about the 
same multiplicity distribution. This has been already 
observed in thermal models but never understood in 
physical terms. In the calculation presented here we 
have used e = Q.15GeV. 

Comparing our completely microcanonical results 
with that of the canonical models |@, 0, we see that 
both give about the same mass yield distribution. 
The strangeness suppression factor which is used in 
|fl"o| is a simple parameterization for the fact that 
strange particles are heavier than their non strange 
counterparts. It depends on the mass difference and 
hence is not identical for all hadrons. This is in agree- 
ment with data. 

In view of these results we can only draw the con- 
clusion that at high energies (i.e. if the average string 
energy is considerably larger than the masses of the 
hadrons considered, as at RHIC energies) the func- 
tional form of the multiplicity as a function of the 
hadron mass cannot serve as a "smoking gun" for 
the plasma. A e + e~reaction, in which certainly no 
plasma but a string is created, produces the same 
functional form. Hence the multiplicity distribution 
cannot tell whether the strings, which are created in 
the AA reactions as well, have interacted to form a 
QGP or whether they will disintegrate independently 
without having it formed. Thus the quest for the 
QGP is more challenging than thought. 
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